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^ . I investigate the discrepancy between the evolution and pulsation masses for 

""1. Cepheid variables. A number of recent works have proposed that non-canonical 

p |. mass-loss can account for the mass discrepancy. This mass-loss would be such 

Q . that a SMq star loses approximately 20% of its mass by arriving at the Cepheid 

. instability strip; a IAMq star, none. Such findings would pose a serious challenge 

. to our understanding of mass-loss. I revisit these results in light of the Padova 
stellar evolutionary models and find evolutionary masses are (17 ± 5)% greater 

^ , than pulsation masses for Cepheids between 5 < M/Mq < 14. I find that 

CN ! mild internal mixing in the main-sequence progenitor of the Cepheid are able to 

^ '. account for this mass discrepancy. 

^ I Subject headings: Cepheids:pulsation stellar:evolution 

oo : 

O ■ 

• S Cepheid variables are an essential step in our determination of extragalactic distances. 

Apart from their use as distance indicators, the regularity of Cepheid pulsation provides 
a well-defined set of observational parameters with which to probe the course of stellar 
evolution for intermediate mass stars. 

Over the preceeding decades much effort has been devoted to reconciling mass determi- 



Introduction 



nations for Cepheids from the various methods at our disposal (see ICoxl (119801 ) for a review). 
Th e longe s t stan ding of these, the Cepheid pulsation mass discrepancy was first revealed 
by IStobid (119691 ) who showed that pulsation masses for bump Cepheid^j were significantly 



^Bump Cepheids are distinguished by a secondary local maximum of the l ight curve seen in Cepheids 
with periods in the range of 6 to 16 days |Bono et allboOObl : lHertzsprun3ll926[) . 
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lower than those predicted by stellar evolutionary models. lAndreasenI (119881 ) showed that by 
artifically enhancing the opacity over the temperature range 1.5 x lO^K < T < 8 x lO^K by 
a factor of 2.5 it wa s possible to remove the discrepancy. More detailed modelling of opac- 
ities by the OPAL (IRogers &: Iglesiasl Il992l ) and the Opacity Project ( ISeaton et al.l [l99J) 
confirmed an increase in this temperature range due to metal opacity. 

The implementation o f new opacities largely resolved the bump Cepheid mass discrep- 
ancy (IMoskalik et al.lll992l ). However, despite convergence, a number of subsequent studies 
have sho wn that the discrepancy remains signif icant an d requires explana t ion. Studies of 



2002, 2006 



Galactic (Natale et al 



Bono et al. 



20021) and SMC flKeller fc Wood 



20071: iCaputo et all 120051) . LMC fjWood et al.l 



1997; Keller fc Wood 



20061 ) Cepheids have shown that 



the masses determined via pulsation modelling are ~15-20% less massive than those ex- 
pected from evolutionary models. Dynamical masses for Cepheids ar e difficult to obtain 



given the low spat ial density of Cepheids. The works of [Benedict et al.l (120071 ): lEvans et al. 



( 120071 . 12OO6I . Il998l ) present dynamical masses for six Cepheids. Albeit with large associated 
uncertainties, these results confirm the conclusions drawn for pulsation modelling; that the 
evolutionary masses appear ~ 15% larger. 

From an evolutionary perspective, Cepheids are understood to be post-red giant stars 
crossing the instability strip on so-called blue loops following the initiation of core-He burn- 
ing. To ascribe an evolutionary mass one takes the Cepheid's luminosity and, using a mass- 
luminosity (M-L) relation that is derived from evolutionary models, derives the mass of the 
Cepheid. The M-L relation can be modified substantially by the treatment of internal mixing 
and mass-loss. Both processes feature complex hydrodynamical and radiative mechanisms 
for which we, at present, only possess empirical approximations. 

The treatment of internal mixing modifies the size of the helium core established during 
the star's main-sequence (MS) evolution. Overshoot at the edge of the convective core of the 
Cepheid progenitor mixes additional hydrogen into the core and hence increases the helium 
core mass. As a consequence the post-MS evolution occurs at a higher luminosity. Mass 
loss, in an ad-hoc manner at least, offers a mechanism to modify the M-L relation by directly 
reducing the mass of a Cepheid. 

The properties of pulsat ion, on the othe r hand, are dependent on the structure of the 
atmosphere of the Cepheid. iKeller fc Wood! (120061 ) show that the morphology of a bump 
Cepheid light curve is highly sensitive to the mass, luminosity, effective temperature and 
metallicity. Hence modelling of the light curve can be used to determine a pulsation Cepheid 
mass that is entirely independent of stellar evolution calculations. 



In the work of iBono et al.l (120021 ) and ICaputo et al.l (120051 ) it is proposed that mass- 
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loss can account fo r the mass discrepancy between pulsation and evolutionary masses. 



Caputo et al.l (120051 ) also conclude that models that incorporate additional internal mix- 



ing in the vicinity of the convective core are not able to explain the mass discrepancy. In 
this paper we revisit these conclusions and present a scenario for resolution of the mass 
discrepancy. 



2. The Cepheid Pulsation — Evolutionary Mass Discrepancy 



A useful way of expressing the mass discrepancy is to form the quantity AM /Me in 
which AM is the difference between the pulsation mass, Mp, and the canonical (i.e. not 
including CCO) evolutionary mass, Me-, as shown in Figure [H In Figured] we also show the 
effects of the inclusion of mild (Ac=0.5) and moderate (Ac=1.0) CCO. The effect of CCO is 
to raise the luminosity of the Cepheid of a given mass as discussed above. 



Caputo et al.l (120051 ) utilise the BVIJK absolute magnitudes derived bv IStorm et al. 



(2004) using distances from the near- infrared surface brightness (IRSB) method. ICaputo et al 



(120051 ) have used the mass dependent period-lu minosity-color relat ion to determine the pulsa- 
tional Cepheid mass using the M-L relation of iBono et al.l (l2000al ). They find that /\M/Me 



ranges from 20% at M ~ 5Mq to ~0% at M ~ IAMq (see Figured]). Such a trend can not be 
explained by a uniform level of core overshoot. Monte-Carlo simulation using the individual 
quoted uncertainties of the data shown in Figure d] reveals a gradient of (3.4 ± 0.5)%Mq^. 



A caveat to IRSB analysis that underlies the results of lStorm et al.l (12004] ) and lCaputo et al 



(120051 ) is the necessary introduction of the poorly unders tood projection factor, p, that em- 
bodies the effects of limb-darkening ( Nardetto et al. 2006h. and is dependent on the pulsation 
velocity and the species under study (INardetto et al.ll2007l ). p is approximated as a function 
of period, and it has been speculated th at this could introdu ce a period dependency in the 
derived distances. However, the work of iFouque et al.l (12007] ) shows that direct HST paral- 
laxes for a sample of Cepheids agree within uncertainties with the IRSB measures albeit for 
only the five stars with distance determinations from both techniques. 



To explain the mass discrepancy found at lower masses, ICaputo et al.l (120051 ) . propose 
mass-loss from the Cepheid progenitor before, or during the central He-burning phase. The 
implied total mass- loss declines from ~ 20% at M ~ ^Mq and vanishes by M ~ 14M0. Such 
mass-loss is seemingly at odds with empirical estimates of mass-loss rates which show that 



mass-loss increases v^ith. stellar luminosity and radius (lReimerslll975l : Ide Jager et al.lll988 



Schroeder fc Cuntzll2007l ). Furthermore, Caputo et al. conclude that CCO can not account 
for this trend in AM/Mg; with mass since this would lead to unphysical AM/M^^a < 
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(their Figure 9; where /Me,k is the mass at a given Ac) for higher mass Cepheids. 



The Cepheid M-L relation implemented in analysis is critical. In their study. ICaputo et al 



(120051 ) have chosen to utilise a linear M-L relation derived from iBono et al.l (l2000a( ) evolu- 
tionary models (Z=0.02, y=0.28 and Ac=0) that do not incorporate mass loss. This M-L 
relation is shown as the dash ed line in Figure [2l However, as shown by the evolutionary 
models of iBono et al.l (j2000al ) there are significant departures from this linear relationship 
for M > 9Mq (see the dotted line in Figure [2]). 

In the analysis to follow we incorporate the non-linear nature of the Cepheid M-L into 
our analysis. The models of (IBono et al.l l2000al . Z=0.02, y=0.27 and Ac = 0) extend to 
M=12Mq, however the derived pulsation ma s ses re ach to 14.9M0. Let us consider the 
stellar evolutionary sequences of iBressan et al.l (119931 ) that extend to highe r masses. These 
mod els (Z=0.02, y=0.28 and Ac=0) implement mass-loss according to the iReimerj (119751 ) 
and Ide Jager et al.l (119881 ) formulation and are shown in Figure [2] as the solid line. For 
masses less than 8Mp there is clearly a different gra dient between thelBono et al.l (l2000al ) and 
Bressan et al.l (119931 ) M-L, a feature discussed by iBeaulieu et al.l (120011 ). What is common 
to both models is a marked departure from a linear M-L in the mass range 9-1OM0. This 
departure results in evolutionary masses that are greater for a given luminosity compared 
with those derived from the linear M-L for masses greater than 9-1OM0. 



The use of the iBressan et al.l ( 119931 ) M-L results in Figure [31 Figure [3] reveals no sig- 
nificant trend in AM /Me as a function of Me (Monte-Carlo simulation of the data and 
associated uncertainties shown in Figure [3] reveals a gradient of (0.3 ± 0.5)%Mq^). This 
demonstrates that a consideration of the non-linear Cepheid M-L accounts for the decrease in 
AM /Me seen in the work of Caputo et al. The absence of a significant gradient also negates 
the argument by ICaputo et al.l (120051 ) that CCO can not be the source of the Cepheid mass 
discrepancy. 



Reanalysis of the results of ICaputo et al.l (120051 ) shows that the evolutionary Cepheid 
mass is (17 ± 5)% greater than that predicted from pulsation modelling. In the next section 
I discuss the merits of a range of possible causes. 



3. The Source of the Cepheid Pulsation — Evolutionary Mass Discrepancy 

The hiding place of the source of discrepancy between pulsation and evolutionary 
Cepheid masses has shrunk dramatically over the last three decades. To account for the 
discrepancy we have three key options. Firstly, it is possible that some input physics in pul- 
sation calculations are not sufficiently described. In this regard, the only input that could 
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affect pulsation to tliis magnitude is tlie description of radiative opacity. Second, tliat tlie 
Cepfieid mass is smaller than their main-sequence progenitors due to mass-loss. And finally, 
that evolutionary calculations underestimate the mass of the He core in intermediate mass 
stars and so underestimate their luminosity. We now discuss each of these possible sources 
in detail. 



3.1. Radiative Opacity 



The pulsation properties of Cepheids are critically dependent on the so-called Z-bump 
opacity arising from the dense spectrum of transitions or iginating from highly io nized Fe. 
The inclusion of th ese transitions in the works of OPAL (IRogers fc Iglesiaslll992l ) and OP 
(ISeaton et al.lll994r) resulted in a substantial increase in opacity at logT fa 5.2. The Opacity 
Project (iBadnell et al.ll2005l ) has included further details of atomic structure (in particular, 
the treatment of atomic inner shell processes) in their calculation of opacity. The new 
opacities do show an increase over the 1992 OP a nd OPAL values of opacity in the Z-bump, 
but at a level of only 5-10% (IBadnell et al.ll2005l ). To account for the mass discrepancy the 
opacity would nee d to be raised by 40 -50%, equivalent to the increase between the early Los 
Alamos opacities (ICox &: Taborlll976l ) and OP and OPAL opacities. Hence the uncertainty 
in radiative opacity is an unlikely resolution to the Cepheid mass discrepancy. 



3.2. Mass-Loss 



The studies of iBono et al.l (|2002| . |2006| ) and lCaputo et al.l (120051 ) propose that mass-loss 
can account for the mass discrepancy. Candidate mass-loss phases include the red giant 
branch phase, sub sequent blue loop evolution, or possibly from the action of pulsation itself 
(IBono et al.ll2006l ). The removal of mass from the Cepheid is a straightforward, albeit ad- 
hoc, way to bring the evolutionary mass in line with that derived from pulsation. The timing 
of, and the changes in stellar structure brought about by significant mass-loss are important 
to the net change in the Cepheid M-L. For instance, significant mass-loss on the MS causes a 
reduction in overall ma ss and hence h elium core mass, resulting in a reduction in luminosity 
in the instability strip (jde Loorelll988l ). Enhanced mass-loss on the red giant branch reduces 
envelop mass, and the material available to the hydrogen-burning shell within t he Cepheid, 



again leading to a reduction in luminosity relative to a star without mass-loss (lYong et al. 
2000l ) . One of the difficulties with the proposal for mass-loss to solve the Cepheid mass 



discrepancy is that it would require the rather artificial bulk removal of material without 
consequent modifications to stellar structure and energy production. 
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Furthermore, standard mass-loss can account for at most a few percent reduction in 
Cepheid m ass and not the 15%-20% required. Mass loss is usually treated the semi-empirical 
relation of iReimerd (Il975l . 'Reimers' law'). While not providing any physical reasoning on 
why the mass-loss is generated, 'Reimers' law' provide s an adequate match to observed 
mass- loss rates over a broad range of stellar parameters (jSchroeder fc Cunta 120071 ). 



Major mass- loss i s exp ected during the red giant branch (RGB) evolution. The models 
of iGirardi et al.l (120001) andlBressan et al.l (119931 ) use a parameterised, empirical fit dM/dt = 
—4 X 10~^^r]L/gR ( Reimera 1975 ). The value of r] is set by a consideration of the masses 
of stars on the horizontal branch (HE) of globular clusters. Determination of rj is made 
difficult due to the variety of HB morphology exhibited by globular clusters. If we consider 
the distribution of effective temperatures for HB stars is entirely due to variable mass-loss 
then 7] must range from to somewhat more than 0.4. Using the canonical value of r/ = 0.4, 
a 5Mq star looses ~O.O3M0 during the RGB phase. To accomodate the mass discrepancy 
seen in Figure [1] this would have to be increased to O.SMq corresponding to a 20-30 fold 
increase in t] which is not plausible. Therefore mass-loss on the RGB does not resolve the 
mass discrepancy. 



Caputo et al.l (120051 ) and iBono et al.l (120061 ) suggest that pulsation may give rise to 



an enhancement of mass-loss. Attempts to measure the mass-los s rates for Cepheids have 



( lDeasylll988l ) and in the radio (IWelch fc Duridll988l ). iDeasy 



been niade u sing IRAS infrare d excesses (IMcAlarv fc Welch 119861) and lUE UV line profiles 



the m ajority of Cepheids of the order of a few x 10 ^ MqJt 



tudy of 


Welch & Duric 


Welch & Durid ( 


19881) 



and iDeasyl (119881 ) conclude that mass- loss during the Cepheid phase is insufficent to explain 
the observed mass discrepancy. In the case of a 5Mq Cepheid, mass loss can account for 
the mass d i screpan cy if either the lifetime in the IS is ten times lon ger than der ived by 



Bono et al.l (j2000al ) or mass loss is thirty times greater than found by iDeasyl (119881 ). Mass 
loss as an explanation for the mass discrepancy in a I2M0 Cepheid is more challenging. The 
mass loss found by Deasy over the entire lifetime of the star can not account for the mass 
discrepancy and a 600-fold increase in the mass loss rate would be required. 



Recently, iMerand et al.l (120071 ) found using near- infrared interferometry that, from a 
sample of four Cepheids all show the presence of some circumstellar material, a Persei, a non- 
variable supergiant residing in the instability strip, does not show evidence for circumstellar 
material, suggesting that pulsation does have a role in enhancing mass-loss. The conversion 
from historical and ongoing mass loss that lead to the circumste ll ar ma terial and the rate 
of mass loss is beyond the scope of the study of iMerand et al.l (120071 ) however. On the 
current evidence I conclude that mass loss does not present a solution to the Cepheid mass 
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discrepancy. 



3.3. Core He Mass and Internal Mixing 

Cepheid luminosity is critically dependent on the He core mass. The mass of the He core 
is determined by the extent of the convective core during core H burning. Classical models 
define the limit to convection via the Schwarzschild criterion. This places the boundary to 
convection at the radius at which the buoyant force acting on a hot clump of material rising 
from the convective core drops to zero. However, the temperature and density regime in the 
vicinity of the convective boundary of the main sequence Cepheid progenitor are such that 
restorative forces in the region formally stable to conve ction are rn i ld, and some signific ant 
level of overshoot of the classical boundary is expected (jZahnlll99ll : iDeng fc Xiongl 120071 ). 



The description of convection is the weakest point in our understanding of the physics 
of intermediate mass - massive stars. Numerical modeling of core convection requires a 
description of the turbulence field at all scales. Three-dimensional hydrodynamical calcu- 
lations capabl e of adequate resolutio n have only recently become a possibility and ar e in 
their infancy jMeakin fc ArnettI 120071 : IPearborn etaPbood : lEgdeton et aPboOsl . I2OO7I 1 In 
the absence of a general theory of CCO, which would enable us to calculate the amount of 
core overshoot for a star of a given mass and chemical composition, a semiempirical phe- 
nomenological approach must be used in calculations of stellar evolution. Several computa- 
tional schemes of various degrees of sophistication in treating the physics of overshoot have 



been discussed in the hterature (iMaeder fc Mevnetl Il988l : iBertelli et al.l Il990l : ICirardi et al. 



2OO0I : IPemarque et al.ll2007l : Istraka et al.lboosh . The most common parameterization of core 



overshoot utilizes the mixing-length formulation, where gas packets progress a distance of 
pressure scale heights into the classically stable region. Ac offers a convenient way to 
parametrize the extension of the convective core it does not constrain its physical origin. In 
addition to the CCO mechanism outlined above, rotationally induced mix ing can similarly 
be inv oked to bring about a si n iilar r ange of internal mixing. As shown by iHeger fc Langei 
(120001 ) and iMeynet fc Maederl (120001 ) mixing in the sheer layer formed at the interface be- 
tween the convective and radiative regions can lead to larger He core masses for massive 
stars. 

At present we must rely on observation for constraint of A^. Observational determi- 



^We quantify core overshoot using the formahsm of iBressan et all (jl98ll ). Note t hat A,, is a factor of 
2.0 times the oversho ot parameter dover/Hp in the formahsm of the Geneva group ISchaller et al.lll99i 
Demargue et aLll2004l ) 
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nations have focussed on populations of intermediate-mass stars (M =5-12M(Tt) , where the 
signat ure o f CCO is exp e cted to be most clearly seen. The studies of iMermilliod fc Maeder 
( 119861 ) and IChiosi et al.l (Il992l ) of Galactic open clusters converge on the necessity for mild 
core overshoot; ft^ 0.5. A number of subsequent stu dies have presented evidence both for 
( iBarmina et al.l |2002| ) and against (ITesta et al.l 119991 ) co re overshoot . From studies of the 
young populous cluster syst ems of the Magellani c Clouds, iKeller et al.l (120011 ) found evidence 



for Ac = 0.62 ± 0.11, while iCordier et al. 



(120021 ) examined the Magellanic Cloud field pop- 
ulation and found the necessity for a level of overshoot between 0.2 and 0.8 pressure scale 
heights. Measurements with the lowest as sociated uncerta i nties are de ri ved from th e pulsa - 
tion modelling o f Cep heid light curves by iKeller fc WoodI kooi . bood ): IBouo et al.l J2002h . 
Keller fc Wood (j2006 ) finds evidence for a weak dependence of Ac with metallicity with Ac 
rising from 0.688±0.009 in the LMC (metaUicity 2/5th Solar) to 0.746±0.009 in the SMC 
(metallicity l/5th Solar). 

While radiative opacity and mass-loss have proven insufficient, increased internal mixing 
remains as the most likely cause of the Cepheid mass discrepancy. The excess of (17±0.05)% 
in evolutionary mass compared to the pulsation mass equates to a uniform level of CCO of 
Ac = 0.67 ± 0.17. This degree of overshoot is within the range of previous studies discussed 
above. The scenario I have outlined offers a straightforward explanation for the findings of 
Caputo et al.l (120051 ) . one that uses canonical mass-loss and mild convective core overshoot. 



4. Conclusions 



In this paper I revisit the conclusions of ICaputo et al.l (120051 ) regarding th e cause of the 



observ ed discrepancy between the evolution and pulsation masses for Cepheids. ICaputo et al. 



( 120051 ) find that Cepheids of 5Mq have 20% less mass, as determined from pulsation analysis, 
than expected from evolutionary calculations, while the evolution and pulsation masses for 
Cepheids of 14M0 agree. In order to explain this finding Caputo et al. propose a scenario 
of non-standard mass-loss:- one that sees increased mass-loss at lower masses and drops to 
negligible mass-loss at M > 14M0. This scheme of mass-loss would be counter to the obser- 
vationally grounded evidence that accumulated mass-loss by the epoch of core-He burning 
increases with increasing stellar mass. Furthermore, Caputo et al. claim that convective 
core overshoot is unable to provide a solution as it can not account for the trend of mass 
discrepancy with Cepheid mass. 

The findings of Caputo et al. are based on a Cepheid mass-luminosity relationship that 
proves erroneous when extrapolated to higher masses. In this paper I show that including a 
full description of the mass-luminosity relation results in a mass discrepancy in which evolu- 
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tionary masses are (17 ± 5)% greater than pulsation masses. The trend of mass discrepancy 
with Cepheid mass is removed. 

I propose that additional internal mixing, as parametrised in the convective core over- 
shooting paradigm, is the primary mechanism giving rise to the mass discrepancy. The level 
of convective core overshoot so derived is Ac = 0.67±0.17, a value which agrees with previous 
determinations from a range of techniques. 

I thank A. Bressan et al. for providing us with unpublished evolutionary models for 
Ac = 1.0. I would also like to thank Peter Wood for discussions during the preparation of 
this paper and Guiseppe Bono for his comments on a draft of this paper. 
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Fig. 1. — Here we show the mass discrepancy as a function of mass. AM /Me expressei> the 
mass discrepancy and is equivalent to the difference between the pulsation mass, Mp, and 
classical evolutionary mass, Me (i.e. Me does not incorporate convect ive core ove r shoot) , 
normalised by Me- Here the extrapolated mass-luminosity relation of iBono et al.l (j2000al ) 
was used (the dashed hne in Figure [21 see text for details). Note the mass discrepancy 
vanishes at higher m asses. Overlaid are the l ocii of models t hat incorporate mild (A = 0.5; 
Girardi et al.l (l2000l )) and moderate (A = 1.0; iBressanI (1200 ll )) convective core overshoot. 
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Fig. 2. — The Cepheid mass-luminosity relation used bv lCaputo et al.l (12005) cons ists of the 
dashed line (from Caputo et al. Equation 2). The models of iBono et al.l (l2000al. Z=0.02, 



y=0.2 7 and = 0) without mass-loss are shown by the dotted line. The iBressan et al. 
(119931 ) mass-luminosity relationship for Cepheids is shown as the solid line. Note the signif- 
icant departures from the linear relation of relation of Caputo et al. 
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Fig. 3. — As in Figure [T] except here the mass-luminosity relation is due to iBressan ei; al. 
(1l993h . The mass discrepancy shows no significant dependence on mass, but rather a uniform 
offset of (17 ± 5)% (see text for details). 



